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Abstract 

This paper deals with faults detection and localization of 
wound rotor induction machines based on principal 
components analysis method. Both, the localization and the 
detection approaches consist in analyzing detection index 
which is established on the latest principal components. 
Once the faults are detected, the affected state variables are 
localized by the variables reconstruction approach. The 
exponentially weighted moving average filter is applied to 
improve the faults detection quality by reducing the rate of 
false alarms. An accurate analytical modeling of the 
electrical machines is proposed and implemented on the 
Matlab software to obtain the state variables data of both 
healthy and faulted machines. Several simulation results are 
prese nte d a nd analyze d. 
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Introduction 

The necessity for having reliable electric machines is 
more important than ever and the trend continues to 
increase. Now, advances in engineering and materials 
science allow building lighter machines while having a 
considerable lifetime. 

Although researches and improvements have been 
carried out, these machines still remain the most 
potential failures of the stator and the rotor. The faults 
can be resulted by normal wear, poor design, poor 
assembly (misalignment), improper use or 
combination of these different causes. Indeed, for 
many years, faults detection in electrical machines has 
been the subject of reflection and research projects in 
various industrial and academic laboratories. 



Several diagnosis and control methods exist and 
already used for the electrical machines monitoring. In 
this paper, Wound Rotor Induction Machines (WRIM) 
faults detection and localization based on Principal 
Components Analysis (PCA) is proposed. PC A is a 
statistical method used for data or state variables 
measurement of systems in operation to monitor their 
behavior. 

The PCA principle consists in reducing the size of the 
representation space of the system [1]. In fault 
detection approach based on PCA, two methods are 
proposed [2, 3], Hotelling's T 2 statistical method and 
Squared Prediction Error (SPE) indicator. The T 2 
statistical is calculated with the "I" first principal 
components while the SPE indicator achieves 
detection with the residual space. However the two 
methods have limitations in faults detection [3, 4]. In 
case of sensors detection, the T 2 indicator is not very 
efficient because the variations due to the failure may 
be masked by normal variations of the variables in the 
first principal components space. And when the 
considered systems are no longer linear, residues 
having high variance contain the modeling errors 
generated by the PCA. Thus the residue having a low 
variance will have less influence on the SPE quantity 
with respect to the residues having a higher variance, 
so that they correspond to the linear redundancy 
relations or quasi-linear. This sensitivity of the SPE 
indicator to the modeling errors can create many false 
alarms. F. Harkat, G. Mourot, and J. Ragot [2] 
proposed a new method for faults detection and 
localization based on the sums of squares of the last 
principal components. 
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For our case, the method proposed by [2, 4] will be 
used for the WRIM faults detection and localization. 

To improve the faults detection and to reduce the false 
alarms, the Exponentially Weighted Moving Average 
(EWMA) filter is used. 

The first part of this article deals with the 
reconstruction principle of the PCA model followed 
by the WRIM modeling. The second part is focused on 
the faults detection method by the detection index (D/). 
The third part is focused on variables reconstruction 
combined with the fault indicator Di for fault location 
aim. The last part is reserved to the applications of the 
PCA approach on the WRIM. Several simulations 
result built with Matlab software are presented and 
analyzed to show the PCA method sensitivity. 

Pea Method Implementation 

PCA methods 

The PCA method is based on a transformation of space 
representation of simulation data. The new space is 
smaller than that of the original space. This method is 
classified as without model methods [5] and can be 
seen as a full-fledged system identification method [6, 
7]. Each variable to be monitored for the state of the 
WRIM are expressed by different units and scales. For 
that, it is preferable to apply a PCA on a centered and 
reduced measures matrix X (columns of zero means 
and units standard deviations) [8]. The orthogonal 
space defined by PCA is generated by the eigenvalues 
and eigenvectors of the matrix correlation R of X. 
These values are sorted in descending order in a 
diagonal matrix. The eigenvalues analysis of the 
correlation matrix R provides information on the 
number of principal components to be retained "I" for 
the PCA model reconstruction [1]. 

The orthonormal projection matrix P formed by the m 
eigenvectors associated with eigenvalues of the 
correlation matrix R is expressed as: 



P = [p l ,p 2 ,...,p m ] 



(1) 



The diagonal matrix A Q f the correlation matrix R 
generated by the eigenvectors associated with 
eigenvalues A sorted in descending order is done by: 



A = diag(A l ,A 2 ,...,AJ 
With \>A 2 >..>A m 



(2) 



The orthogonal matrix which represents the projection 
of X in the PCA new space is T. Mathematically, the 



PCA decomposes X as follows: 



T = XP 



(3) 



T g $R^* m and N is the number of carried out measures 
of variables to be monitored. 

Determination of the structure of the model 

To obtain the model structure, the components 
number "1" to be retained must be determined. This 
step is very important for PCA construction. 
Component number can be determined by using: 



*100> the W 

m 

V k=l J 

With km 

Where the is an user defined threshold expressed as 
percentage. Now, user should retain only the 
components number "1" which was associated in the 
first term of (4). By reordering the eigenvalues, the 
minimum numbers of components are retained while 
still reaching the minimum variance threshold, [9, 10]. 
The vector of principal components is noted by: 



T = [L,t 2 ,...,t m ] 



(5) 



Since the aim of PCA is to reduce the space dimension, 
the "I" first principal components (/ « m) are the most 
significant and sufficient to explain the variability of a 
process. Therefore, the expression of centered and 
reduced matrix X can be written as follows: 



X = X p +E 



(6) 



The matrix X P is the estimated principal part and the 
matrix E the residual part of X which represents 
information looses due to the X matrix dimension 
reduction. They are expressed as follow: 

x P =ip?; (7) 

z=l 

m 

E=2p,T; (8) 

i=l+l 

T is the transpose of the orthogonal matrix. 
Wrim Analytical Modeling 

Fig.l shows the equivalent electrical circuit of the 
WRIM. Each coil, for both stator and rotor, is 
modelised with a resistance and an inductance 
connected in series configuration (Fig. 2). 
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Vj, Ij and Oj (j : A,B, C for the stator phases and a, b, c, 
for the rotor phases) are respectively the voltages, the 
electrical currents and the magnetic flux of the stator 
and the rotor phases, 6 is the angular position of the 
rotor relative to the stator. 




FIG. 1 EQUIVALENT ELECTRICAL CIRCUIT OF THE WRIM 

ij R> ^ 

_J i — i firm . 



FIG. 2 EQUIVALENT ELECTRICAL CIRCUIT OF THE WRIM 
COILS 

Rj and Lj are the resistances and the own inductances 
of the stator and the rotor phases. We note the 
voltages vector ([Vs], [Vr]), the currents vector ([Is], [Ir]) 
and the flux vector ([&s], [Or]) of respectively the 
stator and the rotor: 
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Taking into account the above assumptions, both 
stator and rotor three phase voltages and currents are 
connected to the total magnetic flux by differential 
equations systems [11]. The stator and rotor voltages 
vectors expressions are given by: 



[v.M*,] I'.h^ 

[*s\ = [h] [h] + [M SR ] [I R ] 
M = [L R ] [I R ] + [M RS ] [I s ] 



(9) 
(10) 

(11) 
(12) 



[Rs] and [Rr] are the resistances matrix, [Ls] and [Lr] 
the own inductances matrix, and [Msr] and [Mrs] the 
mutual inductances matrix between the stator and the 
rotor coils. 



Equations (9) and (10) become: 
[v R ] = [R R ] [/«]+ 



dt dt 
d{[L R ][l R ]} d{[M RS ][l s ]} 



dt 



dt 



(13) 
(14) 



By applying the fundamental principle of dynamics to 
the rotor, the mechanical motion equation is [12]: 



t df J V em r 



dG_ 
dt 



(15) 
(16) 

(17) 



Jt is the total inertia brought to the rotor shaft, O the 
shaft rotational speed, [/]=[/a Ib Ic L h Ic]' the current 
vector, f v the viscous friction torque, Cem the 
electromagnetic torque, Cr the load torque applied to 
the machine, 6 the angular position of the rotor with 
respect to the stator, and [L] the inductance matrix of 
the machine. 

Introducing the cyclic inductances of the stator and the 

rotor t =-t and t =-t (Ls is the own inductance 

^sc 2 s RC 2 

of the each phase of the stator and Lr is the own 
inductance of the each phase of the rotor), the mutual 
inductances between the stator and the rotor coils Msr 
and pole pair number p, the inductance matrix of the 
WRIM car be written as follow: 
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(18) 










M SR f 2 


M SR f, 




M s R fi 


M SR f 3 


M s R fi 












M s R fi 


M srA 


M S R f3 





L RC 







M SR f 3 


M s R fi 


M srA 













[L]- 



With 

/^cosipO) 



f 2 =cos(p9 + ^-) 



(19) 
(20) 



264 



International Journal of Energy Science Vol. 2Iss. 6, December 2012 



f 3 =cos(p0-^j-) 

Differential equations system modelling 



(21) 



In choosing the stator and rotor currents, the shaft 
rotational speed and the angular position of the rotor 
relative to the stator as state variables, the differential 
equations a system modeling the WRIM is given by: 



[X] = [AT\[U]-[B][X]) 
With 

[X] = U A h I c l a h I c & 0]', [A ] = 



(22) 
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This model of the WRIM will be used to simulate both 
healthy and faulted operation case of the stator and 
the rotor. The considered faults are resistances values 
increases of the stator or rotor windings due to a rise 
of their temperatures. The following table presents the 
different parameters of the WRIM: 

TABLE I WRIM PARAMETERS 



Symbol 


Parameter 


Value 


Units 


Lsp 


S tato r princ ipal indue tanc e 


0.397 


H 


Lrp 


Rotor principal inductance 


0.397 


H 


Lsi 


Stator leakage inductance 


9.594 


mH 


Lrl 


Rotor leakage inductance 


9.594 


mH 


Msr 


S tato r-ro to r mutual 
inductance 


0,3953 


H 


P 


Numbe r o f po le pairs 


1 




Jt 


Mo ment o f ine rtia 


0.024 


Kg.rrf 


Rs 


Stator resistance 


2.86 


Q 


Rr 


Rotor resistance 


2.756 


Q 


fv 


Viscous frictioncoemcient 


1,444 


mNm/rad/s 



codes and allows us to obtain several matrix data for 
the PCA applications on WRIM faults detection and 
localization. The WRIM is considered faulted from 
t=2s and coupled to a mechanical load at time 2s. Nine 
state variables (m=9) have been chosen to be 
monitored and 10000 measures (N=10000) during 4s 
are considered. 

Fig.3 represents the temporal variations of some state 
variables (Stator current, rotor current, shaft rotational 
speed, angular position and electromagnetic torque) 
showing the steady and transient states of faulted 
WRIM. Fig.4 shows the zoom of the same state 
variables variations but only the part during which the 
machine is in loaded. 



Stator phase "A"current [A] 




FIG. 3 STATE VARIABLES VARIATIONS VERS IE TIME OF 
THE FAULTED WRIM LOADED (STEADY AND TRANSIENT 
STATES) 

Stator phase "A"current [A] 



10 


^0 Rotor phase "a" current [A] 

- 
-10 L 



' Shaft rotational speed [rad/s] 



287.9 
237.85 
287.8 



Angular position [rad] 




10.5 
10. 



Electromagnetic torque [NitT 



The model has been implemented on Mat lab in source 



2.5 2.6 2.7 2.8 2.9 3 3.1 3.2 3.3 3.4 3.5 
Time [s] 

FIG. 4 STATE VARIABLES VARIATIONS VERS IE TIME OF 
THE FAULTED WRIM (STEAD YS TATE) 

Considered faults 

The considered faults are on the resistance values 
which increase due to a rise of their temperature. In 
normal operation, a resistance value variation 
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compared to its nominal value (in ambient 
temperature, 25°C) is faulted machine due to machine 
overload or coils fault [10, 13]. The resistance versus 
the temperature is expressed as: 



R = R (l + aAT) 



(23) 



Ro is the resistance value at To = 25°C, a the 
temperature coefficient of the resistance and AT the 
temperature variation. 

Faults Detection Approach 

Residues generation 

For any measures vector x(k) the equation (6) becomes: 



x(k) = x p (k) + e(k) 



(24) 



x (k) and e(k) vectors represent respectively the 
estimations vector and the estimation errors vector. 

The principal components vector t(k) corresponds to 
x(k) is expressed as: 



t(k) = P'x(k) (25) 
t(k) = [tp(k) te(k)] (26) 

t p g 9t"* z and t e e ^* (m_/) are respectively the " / " first 

principal components vector and the " m-l " last 
principal components vector. 

With this expression (26), there is an similarity on the 
residue vector e(k) and the final components vector 

te(k). 

Detection index "Di" calculation 

The fault detection index is based on successive sums 
of squares of the last principal components [2, 4] and is 
defined as follows: 



j=m-i+l 

i=l, 2,..., (m-l) 



(27) 



At time k, systems are malfunctioning sensing if Di is 
greater than a threshold index noted t? • ol is the false 
detection probability according to the Khi-2 law with 
"m" degree of freedom [14]. One can note a strong 
similarity between the detection index SPE and the 
detection index Dl Indeed, Di corresponds to the SPE 
indicator calculated by PCA model with (m-l) 



principal components. Thus, this threshold detection 
can be calculated with an argument similar to that 
exposed in [4, 14]. 

The process is considered in default at time k if: 



D i (k)>rl 
EWMAfilter 



(28) 



To reduce false alarm and to improve the faults 
detection quality, the EWMA filter is applied at time fc, 
and then the "/nth" filtered vector of the last principal 
components can be written as follow [3, 6 and 15]: 



tAk) = {l-r)t.{k-l) + yt.{k) 



(29) 



y is the forgetting factor (0< y<l) in taking as initial 
condition t ej (0) = and can be calculated by [6, 16] : 

r = l-exp(-l/A r ) (30) 
A t is the time step. 

Finally, the detection index is expressed as: 

m 

a 7 w= Z4w (31) 

j=m-i+l 

And the filtered detection threshold of faults is given 
by [3,10]: 



7 



"if, a 



2-7 



(32) 



It should be noted that many research works uses the 
threshold detection for sensor faults, but our case 
concerns faults detection of systems. 

Faults Localisation Approach 

When a fault is detected, it is necessary to localize or 
identify the involved variables. There are several 
methods for faults localization: 

■ residues structuring approach, 

■ partial PCA approach, 

■ Calculation of variables contributions to the 
detection indicator approaches. 

But [45] showed the disadvantage of the methods 
mentioned above. Then, in this paper, faults location 
of WRIM state variables is based on the variables 
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reconstruction combined with the filtered detection 
index. 

The localization of the WRIM affected state variables 
by the combined approach is based on two methods 
combination: 

■ variables reconstruction by PC A, 

■ detection index 

The method consists in eliminating the fault influence 
on Di when the affected variable is reconstructed. 

To localize faults on the indicator, faults directions 
projected in the residual space should not be collinear 
[8]. After the ; variable number reconstruction, fault 
indicator in noted Dj • 

One can also use the EWMA filter to reduce the 
localization false alarms and to improve the 
localization indicator quality. If D. f is the filtered 

detection index of the ; variable number, the 
localization indicator can be obtained by: 

T if,a 

The variable for which the localization indicator Z^is 

less than one is the offending variable. This method 
can be used for the multiple faults localization in 
reconstructing the supposed faulted variables 
simultaneously. 

Simulation Results and Discussion 

To validate the proposed models and the efficiency of 
the chosen approaches, the established models have 
been implemented in Matlab. Nine WRIM state 
variables (stator three phase currents, rotor three 
phase currents, shaft rotational speed, angular 
position and electromagnetic torque) have been 
considered. The matrixes data of the healthy and 
faulted WRIM obtained by the analytical model of the 
machine are introduced in PCA model to show faults 
detection and localization performances. 

For the electrical machines diagnosis, many methods 
are used to detect the presence or absence of faults, 
occurred at t=2s, and to locate the time when it began 
to appear on the machine windings. Two types of 
faults levels are considered in the system (10%, 30%). 
These values correspond to the rise of the stator or 
rotor coils resistance. We can mention the temporal 



representation (Fig. 5, Fig. 6 and Fig. 7) and the signal 
frequency analysis. Although they have demonstrated 
their efficiency, the state variables representations 
between them also show their advantages. They can be 
performed without mathematical transformation (Fig. 
7). 

Also, the electromagnetic torque variations versus the 
shaft rotational speed clearly show the WRIM 
operation zone in the presence of faults (Fig. 7). After 
several simulations, we suggest some of these 
methods to highlight the place of PCA among them. In 
the Fig.5 and Fig.6, the figures clearly show that it is 
difficult to visualize changes in signals and the fault 
appearance time. 

However, by analyzing the residues of the stator 
current by PCA (Fig.8, Fig.9), the fault appearance 
time is located on the two signals. The case of a 
healthy machine that has a zero residue is almost 
coincident with the x-axis. These observations are 
found in the case of the rotor current (Fig. 8). In Fig.5 
and Fig.6, the presences of faults with the conventional 
temporal representation are no more evident than that 
using PCA method (Fig. 8 and Fig. 9). This one shows 
the residue analysis interest on PCA method. Fig.6 and 
Fig.9 expose the real and residue variations of the 
electromagnetic torque versus time. With Fig.6, the 
fault appearance time is not easy to locate. However, 
with PCA method, the variation of residues in the 
electromagnetic torque with and without faults can be 
easily detected. The real (Fig. 6) and residue (Fig. 9) 
variations of the electromagnetic torque versus shaft 
rotational speed of the machine show again that it is 
much more interesting to treat the state variables of 
the machine with PCA method to detect the presence 
or absence of faults on the windings. The difference 
between healthy and faulted operation (Fig. 9) are 
clearer. It is almost not found in the real variation 
representations (Fig. 6). Fig.8 to Fig. 11 highlight the 
major potential benefits of state variables treatment by 
PCA method which easily shows faults detection and 
locate time of fault appearance. 

With PCA method application, all representation (Fig. 
10 and Fig. 11) well shows the differences between 
healthy and faulted WRIM. In the healthy case, 
residues are zero. When faults appear, the residue 
representations have an effective value with an 
absolute value greater than zero. In the Fig.10 and 
Fig.ll, the healthy case is represented by a right line 
placed on the x-axis. Also, in taking into account "\" 
last principal component, Fig.10 shows peak variations 
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at t=2s (measure number 5000), time at which the 
faults are introduced. The peaks are attenuated 
immediately after but the signals are shifted. 

With "2" last principal components, the detection is 
improved because after the first peaks, other peaks (in 
the presence of faults) are greater than those in the 
case of the healthy machine. 

In the case of the faults detection, Fig.12 and Fig.13 
show the variations of the filtered and the no filtered 
detection index of both "1" and "2" last principal 
components versus the measure number of the faulted 
WRIM. The threshold detection is represented on both 
figures. In the part where the machine is coupled to a 
mechanical load, the shapes exceed repeatedly the 
threshold index. This overrun corresponds to the 
presence of faults. The last principal components 
numbers do not have large influences on the curves 
shape for both filtered and no filtered detection index. 

WRIM Healthy & Faulted: Stator current [A] 



important in faults detection process to avoid alarm 
false. 
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FIG. 5 REAL VARIATIONS VERSUS TIME OF THE STATOR 
CURRENT OF THE HEALTHY AND FAULTED WRM 
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FIG. 6 REAL VARIATIONS OF ELECTROMAGNETIC TORQUE 
VERSUS THE SHAFT ROTATIONAL SPEED OF THE WRM 

However, in the case of the no filtered shape, excessive 
values appear. These values show a bad detection of 
the no filtered data compared to those of the filtered 
data. This behaviour can be corresponding to alarm 
false for some cases. Data filtering is therefore 
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FIG. 7 REAL VARIATIONS OF ELECTROMAGNETIC TORQUE 
VERSUS THE SHAFT ROTATIONAL SPEED OF THE WRM 

Residues :Healthy & Faulted WRIM 







He 


althy 


10% 




30% 
























































> 














Ear 


y fault 


ed 















1.99 1.992 1.994 1.996 1.998 2 2.002 2.004 2.006 2.008 

Time [s] 
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FIG. 9 VARIATIONS OF ELECTROMAGNETIC TORQUE 
RESIDUES VERSUS THE SHAFT ROTATIONAL SPEED RESIDUES 
OF THE WRIM 

Fig.14 and Fig.15 represent respectively the no filtered 
and the filtered localization index versus the WRIM 
state variables. The threshold of the localization index 
is represented on both figures. All variables having a 
localization index founding below the threshold 
variation are the affected variables. 
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Detection index of 1 last principal component 



Dectetion index of 2 last principal components 



2000 3000 4000 5000 6000 7000 8000 
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FIG. 10 DETECTION INDEX OF "1"IAST PRINCIPAL 
COMPONENT VARIATIONS VERSUS THE MEASURE NUMBER 
OF THE WRIM 
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FIG. 11 DETECTION INDEX OF "2" LAST PRINCIPAL 
COMPONENTS VARIATIONS VERSUS THE MEASURE NUMBER 
OF THE WRIM 
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FIG. 12 FILTERED AND NO FILTERED DETECTION INDEX OF 
"1" LAST PRINCIPAL COMPONENT VARIATIONS VERSUS 
MEASURE NUMBER OF FAULTED WRIM 

In the case of the no filtered localization index, only 
variable "4"and "6"correponding to the phase "a" and 
phase "c" rotor currents are not affected. In the filtered 
localization index case, stator phase "a" and phase "V 
are not affected by faults. This last better reflects the 
WRIM behaviour in the case of the considered fault 
type. As in the case of the fault detection approach, the 
use of filter is necessary for faults localization. 
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FIG. 13 FILTERED AND NO FILTERED DETECTION INDEX OF 
"2" LAST PRINCIPAL COMPONENT VARIATIONS VERSUS 
MEASURE NUMBER OF FAULTED WRIM 
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FIG. 14 NO FILTERED LOCALIZATION INDEX VARIATION 
VERS IE THE STATE VARIABLES 
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FIG. 15 FILTERED LOCALIZATION INDEX VARIATION VERSUS 
THE STATE VARIABLES 

Conclusion 

PCA method based on residues analysis has been 
established and applied on WRIM diagnosis. In the 
case of temporal variation and without PCA, the 
electromagnetic torque and the shaft rotational speed 
are the more affected by the considered fault type. An 
accurate analytical model of the machine has been 
proposed and simulated to perform the healthy and 
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faulted data for PCA approach need. 

WRIM faults detection and localization approaches 
based on PCA method are proposed. For that, an 
accurate analytical modeling of the WRIM has been 
carried out. The established models are implemented 
in Matlab. Nine state variables of the machine have 
been considered. Simulation results show the 
efficiency of the detection and localization based on 
respectively the detection index and localization index. 
The use of EWMA filter on both detection and 
localization has helped to avoid some false alarm. Also, 
filtered localization index better show the affected 
variables. 
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